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We present a hybrid model of the unitary-evolution-based quantum computation model and the 
measurement-based quantum computation model. In the hybrid model part of a quantum circuit 
is simulated by unitary evolution and the rest by measurements on star graph states, thereby com- 
bining the advantages of the two standard quantum computation models. In the hybrid model, 
a complicated unitary gate under simulation is decomposed in terms of a sequence of single-qubit 
operations, the controlled- Z gates, and multi-qubit rotations around the z-axis. Every single-qubit- 
and the controlled- Z gate are realized by a respective unitary evolution, and every multi-qubit 
rotation is executed by a single measurement on a required star graph state. The classical informa- 
tion processing in our model only needs an information flow vector and propagation matrices. We 
provide the implementation of multi-control gates in the hybrid model. They are very useful for 
implementing Grover's search algorithm, which is studied as an illustrating example. 
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Q I. INTRODUCTION 

(N 

Q_i A computer is a machine which processes data according 
1^ to a set of instructions. Every computer is a composition 
of hardware on which information is processed and software 
by which information is processed. Hardware is the physical 
O^l part of a computer, while software is a collection of computer 
programs (algorithms) designed to perform a required task. 
c" | A quantum computer (QC) emerges when the computation is 
O executed under the framework of quantum mechanics. There 
j 1 , are several models for quantum computation. 
£h The first conceptual model of a QC, the quantum Turing 
^ machine (QTM), was given by Deutsch pQ. It is rather an 
^.(abstract model and useful for investigating the "computa- 
i__itional complexity." The QTM is a quantum version of its 
classical analogue which allows the superposition of differ- 
ent computational paths. Later on, a rather practical model 



of a QC — the unitary-evolution-based quantum computation 



> 

2^ model (UQCM)— was established [IS]. The UQCM is a 
l quantum edition of "the reversible classical circuit model." 
^— H In the step from classical to quantum, the bits are replaced 
by the qubits, and the logic gates are replaced by the quan- 
turn gates (coherent unitary evolution). Unlike the bits, the 
qubits can exist in a superposition of different computational 
t— I states. Unlike the logic gates, the quantum gates are able to 
^ create and destroy a superposition as well as an entangle- 
• t-h ment. 

The computation in UQCM is run by a sequence of uni- 
^ tary gates and represented by its circuit diagram, where the 
connecting wires stand for the logical qubits or bits which 
carry the information, and the information is processed by 
the sequence of quantum gates. In the end, the result of 
the computation is read out by the projective measurements 
on the qubits. In the well-known textbook by Nielsen and 
Chuang [5], many popular algorithms such as Deutsch's algo- 
rithm [6 , Grover's search algorithm [7], and Shor's factoring 
algorithm [5] are narrated in terms of the UQCM. 

The measurement-based quantum computation model 
(MQCM) is another well recognized model of a QC [H03]. 
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Here, a multi-qubit entangled state — known as a cluster state 
[TP] , or more generally a graph state [TT] — is the main ingre- 
dient, and it provides all the entanglement beforehand for the 
subsequent computation. The computation in this model is 
run by a sequence of single-qubit (adaptive) projective mea- 
surements on the graph state. The methods of MQCM en- 
able one to simulate any quantum circuit on a sufficiently 
large two-dimensional graph state by arranging the spatial 
pattern of measurement bases for the graph qubits according 
to the temporal order of quantum gates in the circuit. 

In the MQCM, the measurements on graph qubits are per- 
formed in a certain temporal order for the purpose of run- 
ning the computation deterministically. Furthermore, the 
measurement outcomes are recorded classically and are used 
for setting the measurement bases for the subsequent mea- 
surements [H]. By contrast, in the UQCM there is no such 
temporal order of measurements, but an order in which the 
unitary gates are executed. 

The task of a QC is to "simulate a quantum circuit." Both 
the UQCM and the MQCM are universal: they can simulate 
any quantum circuit. Where the UQCM uses the unitary 
gates, the MQCM uses the measurements for simulating a 
circuit. In this paper our focus will be on constructing a 
"hybrid" model of the UQCM and the MQCM, which we call 
"hybrid quantum computation model" (HQCM), that com- 
bines elements of both the UQCM and the MQCM with the 
aim of exploiting the strengths of both models. Of course, 
the HQCM is universal as well. There are two main ob- 
jectives of the present investigation. The first objective is 
to develop a theoretical understanding of the HQCM, where 
part of a quantum circuit is simulated by unitary gates and 
the rest by measurements on small graph states. 

The second objective is the investigation of the experimen- 
tal optimization. Both the UQCM and the MQCM possess 
their own advantages along with some similarities. For ex- 
ample, where the implementation of an arbitrary single-qubit 
operation in the MQCM costs a chain of five qubits graph 
state [HI [T3] , it can be implemented rather simply by a uni- 
tary evolution in the UQCM. In return, certain multi-qubit 
gates that are complicated in the UQCM can be realized in 
"one shot" in the MQCM. 

Here is an overview of the paper. In Sec. [Tl] we give a brief 
review of the MQCM, it has two main parts. The first part 
contains the methodology for the computation in MQCM 
with an example and the second part is the classical infor- 



mation processing in the MQCM. We use a portion of Sec. [IT] 
in Sec. IIIII Section IIIII is reserved for a detailed discussion 
on the HQCM, it has three main parts. The first part is the 
methodology for the computation in HQCM and the second 
part is the classical information processing in the HQCM. 
The third part is about the implementation of multi-qubit 
controlled rotations with the HQCM, and this we are going 
to use in Sec. |IV| which is about Grover's search algorithm 
within the HQCM. Finally, we conclude our report in Sec. [V] 
with a summary and outlook. Three appendices deal with 
technical details. 



II. A REVIEW OF THE MQCM 
A. Methodology for the computation in MQCM 

We start this section with a short introduction about the 
preparation of graph states [TU1 E], and then we proceed 
towards the methods for the computation in MQCM [T2"I [T3"] . 
We conclude this section with an example. Throughout the 
text we represent the Pauli vector operator a by (X, Y, Z) 
and the identity operator by /. 

Graph states can be realized in many physical sys- 
tems by first preparing all the qubits of graph Q in an 
eigenstate of their respective Pauli operator X. In other 
words, every qubit 'a' of Q is initialized in the state 
(|0) a + (-l) Ka |l)J/\/2, where n a e {0,1}. Then, entan- 
glement between each pair of nearest-neighbor qubits is es- 
tablished by the controlled- Z gate 

cz(a, b) - |0) a (0| ® + |l) a (l| ® Z<". (1) 

Here, the indices a and b stand for the qubits at the lat- 
tice site 'a' and its nearest-neighbor lattice site of the 
graph Q, respectively. A unitary gate of this kind can be 
generated by turning on the (controlled) Ising-type nearest- 
neighbor interaction for an appropriately chosen time period. 
Experimentally, graph states have been generated by using 
controlled collisions between cold atoms in optical lattices 
[16] and by using linear optics [T7H201 l29 l [30] . 

Mathematically, quantum correlations among the qubits of 
a graph are specified by correlation operators K^'s, which 
are given below. The resultant graph state |3>{ K }) g is an 
eigenstate of these operators, and it is completely specified 
by the set of eigenvalue equations 

K^\$ {K} ) g =X^® ( (g) zA |$ {K} ) e 

\b£nbh(a) / 

= Mr* l*w) s (2) 

with the set of eigenvalues {k} — {n a G {0, 1} | a e Q}. Here, 
nbh(a) stands for the set of all nearest-neighbor qubits which 
are entangled (connected) to the qubit 'a' by the CZ op- 
erations. For every qubit 'a' of the graph state |${ K }) e , 

there exists a correlation operator and an eigenvalue 

K a e {0, 1}. The physical meaning of Eq. ^ is this: there 
exists either a correlation (n a = 0) or an anti-correlation 
(k q = 1) between the outcome of the measurement on qubit 
'a' in the X eigenbasis and the outcomes of the measure- 
ments on all the qubits of nbh(a) in the Z eigenbasis. These 
quantum correlations provide the framework for the compu- 
tation in MQCM. 

Once the resource graph state is ready, then the logical 
qubits — holding the input information — are attached to the 



resource via the same entangling operations given by Eq. ([I]) . 
Unlike the quantum error-correction, here, one logical qubit 
stands for one physical qubit. Now, the computation is car- 
ried out by a sequence of single-qubit (adaptive) projective 
measurements in a certain direction of the Bloch sphere and 
in a certain temporal order. The Bloch sphere offers an ade- 
quate geometrical description for the direction of single-qubit 
projective measurement, where the direction of measurement 
is completely specified by the Bloch vector 

r(9,ip) — (sin 6* cos (,9, sin 9 sin p, cos 9) , (3) 

and the corresponding projector is given by 

P P =(I+(-l) m r-a)/2. (4) 

The measurement outcomes m = and m = 1 mean that 
the measured qubit is projected onto the states with the kets 

It (9, <p)) = cos(#/2) |0) + sin(0/2) 1 1) (5) 

and 

|| (9, <p)) = - sin(0/2) |0) + cos(0/2) |1) , (6) 

respectively. In other words, the choice of measurement basis 
is characterized by the direction of measurement (9, p) in the 
Bloch sphere. 

There are three kinds of measurements in the MQCM 
p"2"l [13] . Measurements along the z-axis effectively detach 
the measured (redundant) qubits from the graph state. Mea- 
surements along the Bloch vector r xy (p) — (cos ip, sin p, 0) — 
it lies in x,y plane of the Bloch sphere — process the infor- 
mation as well as teleport it from one place to another on 
the graph. The significance of this kind of measurements 
is revealed by the example given in Appendix [X] Measure- 
ments along the Bloch vector r zy (9) = (0, sin 9, cos 9)- it lies 
in z,y plane of the Bloch sphere — only process the informa- 
tion. This remark will be illustrated by the example given 
below. 

The two measurement outcomes m = 0, 1 for every qubit 
of the graph state are equally probable because the reduced 
density matrix for each qubit is the completely mixed state 
1/2. In the process of getting the desired operations on 
the logical qubits, one also gets some additional operations. 
These additional operations are called "byproduct opera- 
tors," and they belong to the Pauli group. These byprod- 
uct operators depend on the random measurement outcomes 
and the eigenvalues of the graph state | ) g - The mea- 
surement outcome m G {0, 1} for every graph qubit and the 
eigenvalues {k} are binary numbers, so one can record them 
classically in order to take care of the byproduct operators. 
The classical information processing of this data makes the 
computation deterministic and helps to set the measurement 
bases for the subsequent measurements. Section |IIB| con- 
tains a comprehensive discussion about this matter. Right 
now we are content with illustrating the procedure of MQCM 
with an example. 

Example: The unitary operation for the n-qubit rotation 
around the z-axis is 

U%;;m = exp (-i9Z® n /2) , (7) 

where the superscripts 12. ..n symbolize the logical qubits on 
which this operation will be carried out |15] . One can accom- 
plish this operation by performing a single measurement on 
a (1 + n)-qubits star graph state, the associated star graph is 
shown in Fig. [T|i). In Fig. [lji), the input quantum register 
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FIG. 1: (Color online) (i) This graph is called star graph because 
of its appearance, and the associated graph state |(/>}(i +n ) is given 
by Eq. (jsj. The (blue) circles represent the logical qubits which 
carry the input information; the bonds represent the cz opera- 
tions given by Eq. |TJ, and the ancilla qubit 'o' is represented by 
the (black) diamond. Plot (ii) represents the effect on the input 
state |i/Ji n (n)), when the qubit 'a' of the graph state |<^}( 1+n ) is 
measured in an appropriately chosen basis. 



of n qubits is displayed by the circles and the ancilla qubit 
'a' by the diamond. The input register is in a n-qubit input 
state |"0i n (^)), and the ancilla qubit is prepared in the state 
(|0) a + (^l) Ka |l) a ) /\/2- Then we perform n cz operations, 
represented by the bonds in the figure and given by Eq. ([!]), 
between the qubit 'a' and every logical qubit. In principle, 
all the CZ operations can be performed in "one shot," be- 
cause they all commute with each other. This series of steps 
leads us to the resultant star graph state 



)(l+n) 



V2 



[|0) o ® HM«)) + 

(-l)"-|l) ®(Z«"|^(n)))]. (8) 

The subscript 1 + n reveals that the final graph state is of 
one ancilla qubit and n logical qubits. 

A measurement on the ancilla qubit l a' in the basis 
{| f) I (0, (— l) Ka7r /2)) a } transforms the input state of the 
quantum register into the output state 



|Vout(n)) = (Z°") m ' UltWMnin)). 



(9) 



Here, the direction of measurement lies in z,y plane of the 
Bloch sphere, and m a € {0, 1} is the measurement outcome. 
(iv®™)" 1 " is the byproduct operator, which is represented 
by the dotted-boxes on all the logical qubits in Fig. [ll 11 )- 
After the measurement, all bonds (illustrated in Fig.TlTi)) 
are broken and the qubit 'a' gets projected either onto the 
state | t (6, (~l) Ka7r / 2 ))a (if m a = 0) or onto the state 
| I (6»,(-l) K »7r/2)) a (if m a = 1). We shall use this kind 
of multi-qubit rotations for the HQCM in Sec. |III| 

In contrast to the example R z (tp) given in Appendix [A] 
where the qubits used for input and output are different, in 
the case of U zz ;; z l (6) the input and output states reside in 
the same n logical qubits. In other words, here the informa- 
tion gets processed, but does not get transferred from one 
place to another. As a side remark, the resultant bypr oduc t 
operator (A)"' 1 (Z) Ka in the case of R z (<p) (see Eq. (A3|) 
and the measurement basis {| f,4, (9, (— l) Ka 7r/2)) a } in the 
case of Ug%"£(6) depend on the eigenvalue K a . 

Every quantum gate from the generating set of the Clifford 
group — the CNOT gate 

CNOT(a, b) = |0) a (0| ® I {b) + |1) Q (1| ® X^ (10) 

(the labels a and b are for the control and target qubits, 
respectively), the Hadamard gate 



the 7r/2-phase gate 

R z (n/2) = exp(-Z7rZ/4) 



(12) 



H = (X + Z)/V2, 



(11) 



— can be executed in a single time step in the MQCM |14) . 
This holds because every measurement in these cases is per- 
formed either in the X eigenbasis or in the Y eigenbasis and 
is not influenced by the result of any other measurement. 
Hence all the measurements can be performed simultane- 
ously. A CNOT gate can be achieved with a 15-qubit graph 
state, and both the Hadamard gate and the 7r/2-phase gate 
can be implemented with a chain of five qubits graph state 
[T2I [T5] . Single-qubit and CNOT gates together constitute a 
universal set of gates, and they are realizable in the MQCM. 
In this sense, the MQCM is also universal like the UQCM. 

In order to simulate a complex unitary gate in the MQCM, 
it is customary to first decompose it efficiently into a se- 
quence of elementary gates from the universal gate set. Then, 
the temporal order of gates is transformed into the spatial 
pattern of measurement bases for the graph qubits. After- 
wards, the measurements are performed in the required or- 
der. 

Up to now, we were dealing with the individual gate simu- 
lations only, where we need not to worry about the byproduct 
operators. But in the next section, our focus shall be on the 
simulation of a sequence of gates, where the study of classical 
information processing and the temporal order of measure- 
ments become necessary. Classical information processing is 
needed for taking care of the byproduct operators. A com- 
prehensive discussion about it is provided in the following 
section. 



B. Classical information processing in the MQCM 

This section serves as a summary of the results which were 
discussed in Ref. [14] . When a sequence of gates is simulated 
in the MQCM, the byproduct operator which originates from 
the implementation of gates "passes through" the sequence. 
The propagation of the byproduct operator either transforms 
the next gates in the sequence or the byproduct operator 
in itself gets transformed. The first part of this section is 
about propagation relations for some elementary gates, and 
in the second part we shall define an information flow vector. 
The third part which concludes this section is reserved for 
the propagation matrices for some elementary gate based on 
their propagation relations. We shall use a portion of this 
section for the HQCM in Sec. [Til] 

The structure of the byproduct operator on the logical 
qubit j g {l,...,n} is (IW) 1 '^))'', where Xj and Zj are 
non-negative integers. Both Xj and zj depend on the out- 
comes of measured qubits and the eigenvalues {k} |14j . Their 
dependence on {k} is in our control. For example, the {k} 
dependent parts disappear from the calculation by preparing 
a graph state with k — for all the graph qubits. But we 
cannot control the dependence of the byproduct operators on 
the measurement outcomes which are intrinsically random. 

In Ref. [2], the authors took Xj, Zj € {0, 1}, but we find 
it simpler to take both Xj and Zj as non-negative integers. 
This is permissible because in (X^) Xj and (Z^) Zj only the 
modulo-2 values of Xj and Zj matter. Throughout the paper, 
we reserve the sign '+' for the ordinary addition and the sign 
'©' for the modulo-2 addition. 

In principle, we can correct the byproduct operators — step 
by step — after completing each gate of the sequence under 
simulation. But it is more convenient to choose not to correct 
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them and let them pass through the gates, and just keep 
track of the measurement outcomes in a systematic way using 
simple classical information processing. At the end of the 
computation, either we set the measurement bases for the 
final readout depending on the history of outcomes or we 
just perform the final measurements in the computational 
basis and interpret the result with the help of the record of 
measurement outcomes. 

The propagation of the byproduct operator through a gate 
is given by the propagation relation. The Euler decomposi- 
tion of an arbitrary single-qubit rotation R is 



(13) 



and the propagation relation for R(a, /?, 7) is 

i?(a,/3, 7 ) {Xf{ZY = 

(Xf(^)^((-lfa,(-l)^,(-l) x 7), 



formed to R ((— l) x a, (— l) z f3, (— l) 2 ^), but the byproduct 
operator stays as it is. We can take Eq. (14 1 as an illus- 



tration of the importance of "the temporal order of the mea- 
surements." This is because, when R(a, f3, 7) is a part of a 
circuit, the superscripts x and z are functions of the earlier 
measurement outcomes. In order to determine the right sign 
for the measurement angles a, j3, and 7, we have to wait 
until the necessary measurements are completed [12j. 

Equation (14) also justifies the following points. The 



measurement directions for these qubits lie in x,y plane of 
the Bloch sphere, r = (cos ip, sin ip, 0) with cp i {0,±f}, 
their measurement bases depend on the results of previous 
measurements (see Appendix [A| and Ref. [HI Q3]). When 
ip G {0,±|}, then the directions for +<p and — <p coincide 
and do not get influenced by the outcomes of other mea- 
surements. Measurements of this kind are either in the X 
(ip = 0) or the Y (ip = eigenbasis. The gates from the 
generating set of the Clifford group (given by Eqs. ( 10 )-( 12 1) 
are realized by such measurements. 



The propagation relation for the gate CNOT(a, b) is 



CNOT(a, b) Uf QT = t/s NOT CNOT(a,6), 



(15) 



where 



ttCNOT 

U B 



(x^y°(z^y°{x^) Xb {zWy\ (ie) 



and 



-CNOT 
U Ft 



/j£-(a)\2 a ^(a)y a +ZbfX^) Xa+Xb (Z^y b . (17) 



In case of the CNOT(a, b) gate, Eq. JT5[), the gate stays as 
it is, but the byproduct operator U B gets transformed to 

~ CNOT 

U B . This is also the case for the other two gates from 
the generating set of the Clifford group. The propagation 
relation for the Hadamard gate H is 



H(xy(zy = (xy{z) x H, 



(18) 



and for the 7r/2-phase gate R z (w/2) it is (up to a phase factor 

±0 



R z (n/2)(xy(zy = (xy(zy +x R z (7r/2) 



(19) 



The propagation relations ([15]), ( |18[ ), and (19) can also be 
understood from the definition of the Clifford group which 
maps the Pauli group into itself under conjugation. 



Now, let us define an information flow vector. At every 
stage of the computation, the accumulated byproduct op- 
erator U b upon the logical qubits 1, ...,n is of the form 
n" =1 (^ (i) )^ (ZV))*i . After the implementation of a gate, 
only the values {xj} and {zj} get changed, and the new 
values determine the measurement bases for the subsequent 
gates. These values are processed by a classical computer. 
There is a one-to-one correspondence between the byprod- 
uct operator U b (ignoring the global phase ±1) and a 2n- 
component "information flow vector" X, which is given as 
follows: 



u B = l[(x^yj(zW) 



1 = 



where 



(20) 



(21) 



Here, the multiplication of byproduct operators (up to a 
phase factor ±1) corresponds to the component- wise addi- 
tion of information flow vectors. The information flow vector 
X keeps track of the sign(s) of the measurement angle(s) for 
a gate. In accordance with Eq. (14), the signs of the mea- 
surement angles for the operation R^ (a, /3, 7) on the qubit 
j are determined by the current value of Xj and Zj in X. The 



(14) 




(xi\ 




f Zl \ 


trans- 




\X n ) 




\z n ) 



propagation relations (15), (18), and (19) suggest that none 



of the gates from the generating set of the Clifford group gets 
altered under the propagation of the byproduct operator. So, 
the measurement angles for these gates are independent of 
the values stored in X. 

We can also define a 2n x 2n "propagation matrix" C (g) 
for a gate g, representing the transformation in the informa- 
tion flow vector when the corresponding byproduct operator 
passes through the gate g. The propagation matrices given 
below are derived from the propagation relations (14), (15), 
( [18] ), and (19) with the help of the one-to-one correspon- 
dence given by Eq. (20), and the entries in the information 



flow vectors and the propagation matrices are given only for 
relevant qubits. For the case of n logical qubits, the propa- 
gation matrices for the R-, CNOT-, H-, and i? 2 (7r/2)-gate are 
given in Appendix [B| 

The byproduct operator passes through an arbitrary 
single-qubit rotation R(a, f3, 7) without getting transformed. 
So the information flow vector stays as it is, 




(22) 



C(R) 



The information flow vector gets transformed when the as- 
sociated byproduct operator passes through the CNOT(a, b) 
gate in the following way: 




'10 
110 
11 

v0 1 




(23) 



C(cnot) 



Under the one-to-one correspondence given by Eq. ( 20 ) , the 



propagation relation ( 18 ) for the Hadamard gate H becomes 

(24) 
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and the propagation relation (19) for the 7r/2-phase gate 
R z (n/2) becomes 



x 

z + x 



1 
1 1 



C(««(tt/2)) 



(25) 



The information flow vector and the propagation matrices 
provide a simple description, and they are easily handled by 
a classical computer. 

As a side remark, the temporal order of the measurements 
does not typically follow the temporal order of gates in a 
circuit which we want to simulate with the MQCM. In the 
MQCM, there exists an efficient measurement scheme where 
measurements are performed round by round, and in each 
round all the measurements are executed at the same time 
[14j . The information flow vector is updated after every 
round. After the final round, the result of the computation 
is interpreted from the cc-part of the information flow vec- 
tor X x . An extended discussion of this efficient measurement 
scheme is provided in Appendix [Cj 



III. HYBRID QUANTUM COMPUTATION 
A. Methodology for the computation in HQCM 

In this section, we formulate the methodology of the hy- 
brid quantum computation model (HQCM) by combining 
the advantages of the UQCM and the MQCM. Before going 
into the details, let us first focus on what benefits we can 
get from each of the models in different situations. Here, we 
consider the preparation of a graph state, the set of elemen- 
tary gates for the HQCM, and the simulation of a quantum 
circuit with the HQCM one by one. 

The very first experimental step in the MQCM is the 
preparation of a resource graph state, whereas in the UQCM 
no such preparation of a resource is needed. While preparing 
a graph state, in principle, the initialization of every graph 
qubit in the X eigenbasis can be completed in one shot. To 
this end, we have to "talk" to every graph qubit simulta- 
neously. Consequently, this requires a lot of experimental 
resources, and the very many interactions are difficult to con- 
trol. Likewise, the subsequent two-qubit entangling opera- 
tions (cz(a, b)'s given by Eq. Q) to create the resource graph 
state can be performed in one step, because they commute 
with each other. And, the larger the graph state, the more 
difficult it is to prepare and control the state and to protect 
it against decoherence. So, for the HQCM, we choose not to 
prepare the whole two-dimensional universal graph state at 
once, but instead prepare small (non-universal) graph states 
step-by-step as we need them when the computation pro- 
gresses. Only the star graph states, such as \<f))(i+n) given in 
Eq. |[|J), are required for the HQCM. 

We choose single-qubit operations, the CZ gate and the 
multi-qubit rotation around the z-axis '.'.'"($) f° r an arbi- 
trary value of 9 given by Eq. Q as the elementary gates for 
HQCM. In analogy to the procedure for the UQCM, first we 
"efficiently decompose" any big unitary gate we are trying to 
simulate into a sequence of elementary gates in such a way 
that the number of elementary gates grows polynomially with 
the number of logical qubits, and then every elementary gate 
is implemented one after another. Every single-qubit opera- 
tion and the CZ gate are carried out by the unitary evolution 
under the formalism of UQCM. The ro tatio n U™;"$(0) is im- 



plemented by the method given in Sec. II A (see the example) 



under the formalism of MQCM. The motivation behind these 
choices is explained in the following. 

The implementation of an arbitrary single-qubit rotation 
in the MQCM costs us at least a chain of five qubits graph 
state and four measurements [12l [13]. But it can be real- 
ized quite simply by the unitary evolution of the respective 
single qubit. Furthermore, the Euler decomposition for an 
arbitrary single-qubit rotation R(a, /3, 7) given in Eq. (131 is 
not needed. 

The CZ operations themselves are part of the experimen- 
tal setup for constructing the graph states, and for this we 
have to execute them by the unitary evolution. That is 
why we consider the CZ gate as an elementary gate for the 
HQCM. Furthermore, it is more economical to implement 
CNOT(a, b) by the unitary evolution H^cz(a, b)H (b ^ instead 
of first preparing a 15-qubit graph state and then implement 

it with the MQCM [HH3]. 

Although we already have the universal set of gates (single- 
qubit and CZ gates), we prefer to include U^y^(9) as an 
elementary gate in the HQCM. This is because of two rea- 
sons. The first reason is the optimization. The resource 
(1 + n)-qubit graph state \(/>)(i+ n ) (given by Eq. ^) needed 
for the implementation of is relatively easy to cre- 

ate experimentally. It has only one ancilla qubit, and the 
entanglement can be established in one shot. Furthermore, 
a single measurement on the ancilla qubit is enough to real- 
ize Ul1'"z{9) all together on n logical qubits. While it is also 
possible to decompose the rotation t^'"™ (0) in terms of the 
gates from the universal gate set and implement it under the 
formalism of UQCM, its implementation there will not be so 
optimal, and we cannot regard it as a single unit. 

Generally, one is using either unitary evolution (UQCM) 
or measurements on the graph state (MQCM) in order to 
simulate a quantum circuit. So, the second reason for in- 
cluding Ug%"£(6) as an elementary gate in the HQCM is 
to investigate a model of QC in which a part of a circuit 
(Uzz'.'.'.zi^) rotations) is simulated by the measurements and 
the rest by the unitary evolution, and to see how the classical 
information processing plays its role in such a model. The 
importance of looking at the classical information-processing 
parts is explained in the following. 

Now let us consider the simulation of a sequence of gates 
with the HQCM in which the classical information processing 
becomes crucial. The classical processing does not come into 
the picture of UQCM where the measurements are used only 
for the readout of the final result of computation. In all those 
schemes where measurements are needed for the computation 
(for example quantum teleportation [5]), the classical infor- 
mation processing in parallel is essential. In the HQCM also, 
classical information processing is needed, because the rota- 
tions UHy^{9) are executed by the measurements. But here 
the classical information-processing parts are rather simple 
and straightforward, requiring only the information flow vec- 
tor and the propagation matrices. A comprehensive discus- 
sion of this is given in the following section. 



B. Classical information processing in the HQCM 

In this section, our focus will be on the classical 
information-processing parts of the HQCM. We only need 
the information flow vector and the propagation matrices 
for the elementary gates. We first redefine the information 
flow vector in the context of HQCM, and then discuss the 
propagation relations as well as the propagation matrices for 
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the elementary gates. 

At every computation step r, the form of the byproduct 



operator is the same, as given in Eqs. (20) and (CI) 



Ub(t) 



3=1 



(26) 



So, the form of the related information flow vector 1(t) = 
I(xj(r), Zj(r)) is also the same as given by the 2n x 1 column 
vector in Eqs. (pOl and (|2lj). But in the HQCM, there are 



some differences in comparison to the efficient measurement 
scheme of MQCM given in Appendix [C] In that scheme of 
MQCM, the index T of l(t) stands for 'the measurement 
round.' But in the HQCM, every elementary gate is taken 
as a single computational step, and the index V of 1(t) is 
the label for them. In the MQCM, l(t) gets updated after 
each round, but in the HQCM it is updated after each gate. 

In this scheme of MQCM, the initial value of the in- 
formation flow vector 2w? t is determined by the set of 
eigenvalues {n} and by some particular gates. But in the 
HQCM, just before starting the computation all the entries 
of 1(0) = 1?£™ are zeros, i.e., both Xj (0) = and ^(0) = 
for all j — 1, 2, • • • , n. That means that the byproduct op- 
erator at t — is the identity operator / on every logical 
qubit. In fact, the first relevant byproduct operator appears 
in the computation when the first multi-qubit rotation is im- 
plemented, and then the information flow vector gets some 
nonzero entries. In the MQCM, the information flow vector 
gets updated from l(t — 1) to l(t) after the tth measure- 
ment round. In the HQCM, the information flow vector gets 
updated from 1(t — 1) to 1(t) after the implementation of 
rth gate. 1(t) influences the (r + l)th gate of a quantum 
circuit under simulation. Similar to the case of the UQCM, 
the total number of computation steps (the logical depth) is 
denoted by r max , which is the total number of elementary 
gates used for the computation. Furthermore, r max is also 
the total number of steps taken by a classical computer for 
the classical information processing in the HQCM. 

Let us turn to the issue how one can interpret the final re- 
sult of the computation in the HQCM. In the case of UQCM, 
every gate of a circuit is executed by their respective unitary 
evolution, and the final readout measurements are performed 
in the computational basis. In this case, the output state 
|out) gets projected onto the state |Muqcm) 
after the final readout measurements, 



9>=1 \ s l) 



IM 



UQCM/ 



n 



iU) + (-i)^zb) 



lout) , 



(27) 



where Sj € {0, 1} are the readout measurement outcomes for 
the logical qubits j = 1, 2, • • • , n. 

In the case of HQCM, the final state of the output reg- 
ister will be U B(r max ) |out) after performing the last gate 
of the same circuit. Without loss of generality, like above, 
we consider the computational basis for the final readout, 
where Sj £ {0, 1} are the readout measurement outcomes 
for the logical qubits j = 1, 2, • • • ,n. That means that the 
output state J7s(T max ) |out) gets projected onto the state 



Af, 



HQCM/ 



after the readout measurements, 



IM 



HQCM/ 



11 9 t/i 



(r max ) |out) . (28) 



We can transform Eq. ( 28 ) with the help of Eq. ( 26 1 into 



I Mi 



HQCM/ 



) n — i° ut > • 



(29) 



The inference we get by comparing Eq. (27) and Eq. (29) 
is: the readout measurements on the state out) with the re- 
sults {Sj} give the same circuit-output as the readout mea- 
surements on the state C/s(r max ) |out) with the results {sj}, 
and these sets of results are related by 



s 3 = sj + Xj(T max ) for all j G {1, 2, • • ■ , n} . 



(30) 



That is how one can interpret the final result of the compu- 
tation with the help of l x (Tmax) in the HQCM. 

Let us turn to the propagation relations for the elemen- 
tary gates. An arbitrary single-qubit rotation around an axis 
r(0,ip) (defined in Eq. p|) by an angle a is 



Rp(a) — exp i—f.a^j 



(31) 



The byproduct operator passes through this gate without 
any change, but it changes the axis of rotation of the gate 
from r to r'. The propagation relation for R?(a) is given by 



Rp(a)(X) x (Zy = (Xy(Z) z R?,(a), 



(32) 



where 



r ' = ((-I) 2 sin6» cos<p : (-l) x+z sin6» sinv?, (-If cos6») . 

(33) 

In other words, the angles 9, ip that define the axis of rotation 
r get transformed as 6 — > (xit ~ 9) and tp — > (— l) x (zn + if). 
The byproduct operator passes through R?(ot) without get- 
ting transformed, that means that the propagation matrix 
C(R) is the same identity matrix as given in Eqs. (22) and 
(Bl|. Every single-qubit unitary operator in SU(2) follows 
this propagation relation, and it becomes the propagation re- 
lation ( 18 ) for the Hadamard gate when 9 = ip = tt/2 and the 
propagation relation ( 19 ) for the 7r/2-phase gate when 9 = 0. 



So, the Hadamard- and the 7r/2-phase-gate remain special 
cases in the sense that the byproduct operator changes un- 
der the propagation but not these gates. They are executed 
by the unitary evolution like any other single-qubit gate, but 
for the classical information-processing parts of the H QCM 
we shall use their propagation matrices given by Eqs. (B3) 
and (iBl 



The propagation relation for the next elementary gate 
cz(a, b) (defined in Eq. p])) is 



where 



cz(a, b)U C s = U B cz(a,b), 



U c ^ = (lW) I «(2(»)) 2 «(l( t )) I »(2( , '))^ i 



and 



U, 



+x a 



(34) 



(35) 



(36) 



Under the one-to-one correspondence given in Eq. ( 20 ) the 
propagation relation (34) becomes 




'1 0' 
10 
110 

a o o i 

V 

C(cz) 




(37) 
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When the control qubit 'a' and the target qubit belong to 
the set of n logical qubits [a ^ b), then the propagation ma- 
trix C(cz(a, b)) can be generated by the following relations: 



[C ra (cz(a, b))] kl = [C 2Z (cz(a,6))] fc; = 8 kU 
[C xz (cz(a, b))) kl = S ka Si b + S kb 5i a , 
[C zx (cz(a,b))] M = 0. 



(38) 



Note that Eqs. p7|) and (38 1 are different from Eqs. (p3l and 



(B2|. The CZ and CNOT gates are interconvertible by using 
the Hadamard gate, i.e., H^Cz(a,b)H^ = CNOT(a, b) and 
the same is true for their propagation matrices, i.e., 



C(JjW)C(cz(a,6))C(ffW) = C(cnot(o,6)). 
The propagation relation for t/^f.'.'.'J (0) is 

ult.mu B = U B Ult~{{-lT6), 

where U b is the same as given in Eq. p0| and 

n 

x = Xj . 



(39) 



(40) 



(41) 



In this case, the measurement angle 9 gets modified under 
the propagation, but the byproduct operator stays as it is. 
So, the propagation matrix C(f7*|"'™(#)) will be the iden- 
tity matrix, which can be defined in the same way as the 
C(R) in Eq. (Bl). The inferences we get here are the fol- 
lowing: (1) The Hadamard-, the 7r/2-phase- and the cz-gate 
remain unchanged under the propagation, while the byprod- 
uct operator gets altered. (2) The Single- and multi-qubit 
rotations (with nontrivial angles) get transformed, while the 
byproduct operator stays unaltered under the propagation. 

Now we have all the basic tools for the HQCM. We shall 
first describe some important examples, and then proceed to 
the implementation of Grover's search algorithm within the 
HQCM. 



C. Controlled operations with the HQCM 

In this section we are considering n-qubit controlled rota- 
tions around the z-axis, which are defined by 



A 1 - c C/i c + 1) -™(^) = (I® c - |l-l)i... c ® I® {n - C) + |l...l)i... c ® U( c +V- n (6), 



(42) 



where the qubits labeled 1 to c arc the control qubits and the 
qubits labeled c + 1 to n are the target qubits. Only if every 
control qubit is in the state 1 1) , then the (n— c)-qubit rotation 

Uz°^z n (&) operates on the target qubits. We discuss the 
HQCM implementation for three different values of c: c = 
1 (single-control), c = 2 (double-control) and c = 3 (triple- 
control). 



The single- control gate A 1 !/'^ 



-26) 



First, we decompose A 1 J7 Z " ™(— 29) in terms of multi-qubit 
rotations like U\ 2 "^{9) given by Eq. Q. In order to have 
the logical qubit 1 as the control and the rest as the target 
qubits, we express the n-qubit rotation about the z-axis as 



t£:::?(0) = |o>i<o|®E£ 



1 (0) + |i>i<i|®e£;;?(-0). 



Consequently, we get the required decomposition: 



(43) 



We observe that ■■■"(#) is symmetric under permuta- 
tion of the qubits, so we can take any qubit as the control 
and the remaining qubits as targets. We consider a 'multi- 
qubit rotation about the z-axis' as a single unit, and then 
A 1 Uz ' ™(—29) costs only two units of this kind. 

The circuit representation of Eq. ( 43 ) is given in Fig. [2] In 
Fig. [2](i) , the first and the second rectangular boxes depict 
U\ 2 Z ;;™{9) (the first rotation) and U 2 y z {—9) (the second rota- 
tion), respectively. In practice, both of them are realized — by 
using a single ancilla qubit and a single measurement — with 
the methodology given in Sec. II A (see the example). More- 



over, after executing the first rotation we bring back the an- 
cilla qubit into an eigenstate of X and use it for the second 
rotation. The implementation of U^y"(6) and U z y™{—6) 



require (1 + n)-qubit and n-qubit star graph states, respec- 
tively (see Fig. [iji)), where the ancilla qubit is connected 
to the relevant logical qubits. The eigenvalues of the ancilla 
qubit corresponding to the first and the second rotation are 
K\ and «2; and the measurement outcomes are mi and m2, 
respectively. 

The classical information processing for this case has three 
main parts. The first part deals with the change in the mea- 
surement angles due to the byproduct operator Us.in- t/s,in 
appears just before implementing the first rotation and is 
denoted by the dashed vertical line at the input section in 
Fig.[2] When the gate A x £/f;;;™(-20) in itself is a part of a cir- 
cuit under simulation, then due to the implementation of pre- 
vious gates the byproduct Ub, in has emerged prior to the ex- 
ecution of A 1 U z --™(—29). Without loss of generality, we take 
U B , in = N2 =1 (XV))*i (ZW)* the same as given in Eq. |20|. 
Only the x-part of the corresponding information flow vector 
l x ,in influences the measurement bases {| f, 4 (±#> 7r /2)) a } 
for both rotations. According to Eq. ( |4"o"| ), the angles 9 for 
the first and —6 for the second rotation get altered. 



(-l) x 9 
-{-1) X ~ X H 



for 
for 



u 12 - 

^ zz.. 
jj2... 



(44) 



where x is given by Eq. (41) 



The second part deals with the eigenvalues n\ and «2, 
which influence the azimuthal angle 7r/2 of the measurement 
bases in the following way: 



tt/2 -> (-1) Ki tt/2 for Ulf~™{6), 
tt/2 -> (-l) K2 7r/2 for U%;;£(-9). 



(45) 



The third part manages the contribution of measurement 
outcomes mi and m.2 to the byproduct operator Us^n- The 
implementation of both the first and the second rotation 
cause the byproduct operators {Z® n ) mi and (Z^"" 1 ))™ 2 
on the relevant logical qubits. Furthermore, these byproduct 
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FIG. 2: (Color online) (i) The quantum circuit merely represents 
the temporal order of rotations for A 1 ^ ;.'.'"(— 26). Horizontal 
lines stand for n logical qubits. The first (green) rectangular box 
(from the left) symbolizes the n-qubit rotation U and the 
second box symbolizes the (n — l)-qubit rotation Uj";£(—0). Both 
of them are executed under the scheme described in Sec. II A (ii) 
This circuit represents A 1 C/f;;;"(— 20), where qubit 1 is the control 
qubit and the other qubits are targets. The (red) dashed vertical 
lines at the input and the output section stand for the byproduct 
operators 17b, in and Ub, out, respectively. Circuits (i) and (ii) are 
equivalent. 



operators update Us.m to Ub. out- Ub, out is denoted by the 
dashed vertical line at the output section in Fig. [2] Only 
the z-part of the information flow vector X z j n gets changed, 



while the x-pa.it remains as it is, i.e., X x 



— j 



U B ,out = (x {1) ) Xl {Z {1) ) Zl+mi Y\_{x {i) ) Xi (zU ) ) z i +mi+m * . 

(46) 



2. The double-control gate A 12 t/ Z 3 ;;™(40) 

We have to combine two additional units U^ 3 ^ n (— 0) 
(the third rotation) and f/J.'.'.'" (0) (the fourth rotation) 
with A 1 // 2 "'"(— 20) for the purpose of getting the gate 
A 12 L^;;;™(40). In other words, A 12 17 Z 3 ;™(40) with two con- 
trol qubits 1 and 2, is made up of four rotations, and its 
decomposition is given by 

A i2 ut:m = ut:(o)utt.::(-mt.:(-mit ::((>)■ 

(47) 

Figure [3ji) illustrates the temporal ordering of the multi- 
qubit rotations given in Eq. (47) by the rectangular boxes. 

for A 12 Uf^(4:0) is similar to that of 
All the rotations— (7 12 -™(6») (first), r/f-™(-0) 
Uz 3 :: z n (-6) (third), and U^;;;f(6) (fourth)— are 



The treatment 
Sec 



III C I 

(second), 

performed one after another under the scheme given in 
Sec. " 



II A 



(see the example). After initializing the ancilla 
qubit in the X eigenbasis, we prepare a necessary star graph 
state for the first rotation, and then the ancilla qubit is mea- 
sured in the appropriate basis. The measurement outcome 
is recorded, and the ancilla qubit is brought back again into 
an eigenstate of X (recycled) for executing the next rota- 
tion. In this way, we can use the same ancilla qubit for all 
the rotations. K\, K2, k , and K4 are the eigenvalues of the 
ancilla qubit, and mi, ma, m.3, and 777,4 are the measure- 
ment outcomes corresponding to the first, second, third, and 
fourth rotation. As a side remark, one can also choose to 
perform these four rotations at the same time by using four 
different ancilla qubits, but this would require more hardware 
resources. We can take A 12 [/J* ,"'"(40) as a single time step, 



1 -!- 



2 -L 



i l{e)_{-e)-{-e)-{8)\ 





1 

— T 


1 

T 


-{26) 


— -+ 


(40) J- 
1 




i 1 

-1 _L 


1 

j_ 
1 



u' 



(i) 



(ii) 



Us. m U' Bjs 
(iii) 



FIG. 3: (Color online) The horizontal lines stand for n logical 
qubits. (i) Four (green) rectangular boxes (count from the left) 
represent Ull:;V:(6), U 2 z ;;;l\-6), U™;*" (-0) , and Ut.:"(0), re- 
spectively. Each rotation is realized under the scheme described 
in Sec. II A (ii) (from the left) The first (green) rectangular box 
symbolizes the n-qubit operation A 1 U^',;;"(~26), and the second 
one symbolizes the (n— l)-qubit operation A 1 U^'/'"(28). Both of 
them have the qubit 1 as control, (iii) The diagram represents 
A 12 f7|;;;"(40), where the qubits 1 and 2 are the control qubits. In 
(i) and (iii), the (red) dashed vertical lines at the input and out- 
put section stand for the byproduct operators Ub. in and Ub, out, 
respectively. Circuit (ii) merely depicts the intermediate stage of 
circuits (i) and (iii), and they all are mutually equivalent. 



because, in principle, all the four rotations can be executed 
at the same time. 

The classical information processing for this case also has 
three main parts. The first part deals with the modification 
in the measurement angles because of the byproduct operator 
U BM = Ui=i{ x(j) ) Xi ( zU) Y* , which is represented by the 
dashed vertical line at the input section in Figs. [3](i) and 
[3^ iii). Here also, only I x ^ n influences the measurement angle 
±9 for every rotation. 



(-iy 



(-iy 



if- 
■iy- 



X2 I 



for 
for 
for 
for 



uj;::"{d), 

J2...n { _ 6l 

ulh:z(~e), 

r 3 -™(6»), 



U 2 



(48) 



where x is the same given by Eq. (41 1. 

The second part manages the influence of the eigenvalues 
K\, k 2 , K3, and K4 on the azimuthal angle 7r/2 of the mea- 
surement bases {| til (±0, 7r /2))a} in the following way: 



tt/2 -> (-1) Ki tt/2 for 

tt/2 (-1) K2 tt/2 for 

tt/2 -> (-1) K3 tt/2 for 

tt/2 (-1) K4 tt/2 for 



Uzz::.z(o)> 
ut:"(~o) 

Tjl3...n(_c 

ul 



(49) 



ZZ...A 

r3...n 



(0). 



The third part handles the random measurement outcomes 
mi, 77i2, ma, and TO4; which cause the byproduct operators 
(^«) mi , (£»(«-i)) m2 , (z^™- 1 ))™ 3 , and (Z®("" 2 )) m4 , re- 
spectively, on the relevant logical qubits. Furthermore, they 
change UsAn into UB,ont by their contribution. Ub.oxa is 
denoted by the dashed vertical line at the output section in 
Figs.[3]Ji) and[3pii). So, only the z-part of the corresponding 
information flow vector I z ,i n gets changed to I z , ou t, while 



1, 



/ zi + mi + m 3 \ 
z 2 + mi + m 2 
zz + m 

z n + m 



(50) 
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FIG. 4: (Color online) The quantum circuit for implementing the 4-qubit gate A^Z 6 . From the top, three black horizontal lines 
stand for the control qubits 1, 2, and 3, and the next two gray horizontal lines stand for the work qubits 4 and 5, which are prepared 
in the state |+)® 2 . The black horizontal line at the bottom stands for the target qubit 6. Every 3-qubit gate is the special case of 
A 12 U^;y" (49) (for n = 3 and 8 = ±7r/4), and they are realized by the four rotations according to Fig. |3| The Hadamard gate H and 
cz(5, 6) are executed by the unitary evolution. The computation steps (r) are expressed by (red) dashed vertical lines for the classical 
information processing according to Table [I] 



where 



m = mi + 77i2 + TI3 + TO4. 



We emphasize two points. First, A 12 £/ 3 ' ""(40) is sym- 
metric under permutation of the qubits. So, we can take 
any two qubits as controls and the rest of the qubits as tar- 
gets by using only four multi-qubit rotations. If we con- 
tinue along the same direction, then the controlled rotation 
A 1 ~ c Uz°^z'" n (6) with c control qubits requires 2 C units of 
rotation. The number 2 C is independent of the number of 
target qubits (77 — c), but when c becomes order of n, then 
2° becomes exponential in n. In order to fix this exponential 
growth problem, we need some extra work qubits [3 . Our 
next example will justify this remark. 

Second, A 12 C/ Z 3 ; ™(46») becomes A 12 t/ 3 (4(9) for 77 = 3. The 
gate A 12 U^(49) is equivalent to the Deutsch's universal gate 
A 12 [ji? 3 (4#)] [5] up to single-qubit gates, provided the an- 
gle '40' is incommensurate with ir. We are going to use 
A 12 [/ 3 (±7r) with two work qubits for implementing our next 
gate A 123 Z 6 . 



3. The triple- control gate A 123 Z 6 

We put the bits and pieces of the HQCM together as a 
summary by taking the gate A 123 Z 6 as an example. The 
complete scheme about its implementation in terms of its 
circuit diagram is shown in Fig. |4j and the associated classi- 
cal information- processing parts are given in Table [i] 

First, we efficiently decompose the gate A 123 Z 6 into a se- 
quence of elementary gates (single-qubit gates, the CZ gate, 
and C/ 12 '" z (#)). The temporal order of the elementary gates 
for A 123 Z 6 is depicted by the circuit diagram in Fig.|4j where 
the qubits 1,2, and 3 act as the control qubits and the qubit 
6 as the target qubit, and they all are represented by the 
black horizontal lines. The work qubits 4 and 5 are initial- 
ized in the state |+)® 2 , and they are represented by the gray 
horizontal lines in the figure. The work qubits are used to 
make the decomposition of A 123 Z 6 economical. 

The 3-qubit gates A 12 [/ 4 (tt) (1st gate), A 34 [/ z 5 (tt) (3rd 
gate), A 34 J7 2 5 (-7r) (7th gate), and A 12 [/ 4 (-7r) (9th gate) are 
represented by rectangular boxes with double-control. Every 



3-qubit gate is further decomposed into four rotations around 
the z-axis according to Eq. (47), here n = 3 and 9 = ±7r/4. 
Furthermore, each rotation is executed by preparing a re- 
quired star graph state, followed by the measurement in the 
appropriate basis. The detailed methodology is already men- 
tioned in Sec. |III C 2| The Hadamard gates H are displayed 
by rounded rectangles, and the two-qubit gate cz(5,6) is 
shown by the rounded rectangle on the qubit 6 with the 
qubit 5 as control in Fig. |4j The Hadamard and the cz(5, 6) 
gates are executed by the unitary evolution. 

The classical information-processing parts for A 123 Z 6 are 
handled by a classical computer according to Table [I] In this 
table, the first column is for the computational steps t, which 
are represented by the dashed vertical lines in Fig. |4j There 
are ten vertical lines in the figure and ten rows in the ta- 
ble for the ten computational steps from to 9. At each 
vertical line the information flow vector T(r) gets updated. 
The second and the third columns are reserved for X x (t) and 
Z z (t) respectively. If required, the change in the measure- 
ment angles for the next gate based on the updated value of 
I(t) is calculated; they are given in the fourth column. After 
performing the measurements in the appropriate bases, the 
measurement outcomes are recorded in the fifth column. 

Let us go through the table row by row. Before starting 
the computation (in the first row r = 0), all the entries of 
both I x (0) and T Z {Q) are zeros (initialization). So, there is 
no change in the measurement angle for each of the four rota- 
tions associated with the gate A 12 f/ 4 (7r) (1st gate). The mea- 
surement outcomes mn, 77712, 77113, and mu corresponding 
to the four rotations U™(9), U%t(-0), UT£{-9), and U*(9) 
are recorded. These outcomes give some nonzero entries to 
X z (l) according to Eq. (50 1. The measurement outcome rrijk 



corresponds to the feth rotation of the jth 3-qubit gate. The 
next gate in the circuit is the Hadamard gate H^ 4 \ which 
does not change under the propagation of the byproduct op- 
erator; that is why the forth column of the second row r = 1 
is empty. The iJ-gate is realized by the unitary evolution, 
that is why the fifth column of the second row is also empty. 
But the ff-gate changes the information flow vector 1(1) into 
1(2) under the propagation relation given by Eq. (18), and 



the p ropa gation matrix for the H-g&te is given by Eqs. (24) 
and (|B3|. The 3rd gate is A 34 C/ z (7r). The measurement 



angles ±9 only for the rotation J7 ZZZ (0) and J7 ZZ (— 9) get in- 
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TABLE I: The classical information-processing parts for A Z 



Z*(t) 



Angle ±0 (here 6» = tt/4) 



Measurement outcomes 



M 






W 
M 





W 

( o \ 



mi 


V o ) 
( o \ 




mi 


V o / 
/ o \ 

{) 
{) 

mi 

V / 

( o \ 



mi 
m 3 

V / 
/ \ 



mi 


V o / 
( o \ 



mi 


V o / 





m 3 l + "132 + m 7 l + m 72 


V ) 
I \ 




m 3 i + m 32 + m 7 i + ?n 72 


V o / 



J 



J 



M 






\oJ 

/ mil + mi 3 
run + "112 


mi 


V o J 

1 'mil + mi3 ' 
mn + »"12 




V o 

I mn + mis\ 
mn + mia 
m 3 i + m 33 
m 3 l + "»32 
m 3 

V o 

''mil + mi3 ' 
mil + mia 
m 3 i + m 3 3 
m 3i + m 32 


V 
I mn + mi 3 

mn + mia 
m 3 i + m 33 
m 3i + m 32 


V m 3 
f mn + mi 3 ^ 
mn + mi2 
m 3 i + m 33 
m 3 i + m 32 

m 3 

V m 3 
mn + mi 3 
mn + mia 

m 3 l + m 33 + m 7 i + 
m 3 i + m 32 + m 7 i + 
m 3 + m 7 
m 3 
mn + "il3 
mn + mi2 
"i 3 l + m 33 + m 7i + 
mi 
m 3 + m 7 
m 3 

/ r«n + mi3 + mgi + 
r«n + mi2 + mgi + 
m 3 i + m 33 + m 7i + 
mi + m 
m 3 + m 7 
V m 3 



m 73 
m 72 



No change in angle for 

No change in angle for 

No change in angle for 

No change in angle for 



u 2 z i(-o) 

Ult(-9) 
Ui(9) 



-> (-l) m i0 for Ulf z {6) 

-0^-(-l) m i0 for u^l(-e) 

No change in angle for U zz ( — 0) 

No change in angle for U z (9) 



-0 ->■ -(-l) m i0 for Ulf z {-0 

-> (-l) m i0 for Uf> z (0) 

No change in angle for U zz (0) 

No change in angle for U z (—9) 



m 93 
m 92 
m 73 



-e->-(-l)™6 for (7^f(-i 

e -4- for 

-4- (-1)™6> for (7^(9) 

-0-y-(-l)™0 for !7*(-0) 
m = m 3 i + m 32 + m 7 i + m 72 



mn for UH*(0) 
mia for !7jJ(-0) 
™13 for (7^(-e) 
m u for U*(6) 
mi = mn + mi2 + mi 3 + m i4 



m 3 i for !Z^J(0) 

m 32 for Utli-6) 

m 33 for 

m 34 for U r i(0) 
m 3 = m 3 i + m 32 + m 33 + m 34 



m 7 i for U^f z (-0) 

m 72 for C/«(e) 

m 73 for UH(0) 

m 74 for C/f(-0) 
m 7 = m 7 i + m 72 + m 73 + m 74 



mgi for UlHi-0) 

m 92 for ajj(e) 

m 93 for Ult(0) 

m 94 for U z (-0) 
m 9 = m 9 i + m 92 + m 93 + m 04 



fluenced by I x {2) according to Eq. ( |48| ). The measurement 
outcomes mai, 777,32, 777.33, and 777.34 only transform I z (2) into 
I z (3). In this way going through Table [I] along with Fig. [4] 
explains the whole scheme, and the final output result is in- 
terpreted according to Eq. (30) with the help of I x (9). 



Here, the z-part of the information flow vector I z (t) gets 
the new entries from the implementation of 3-qubit gates 
only according to Eq. (50 1. The entries of I(t) get ma- 
nipulated under the propagation of the byproduct operator 
through the Hadamard gates and the cz(5, 6) gate according 
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to the propagation relations (18) and (34), respectively. But 



the propagation of the byproduct operator does not change 
the H- and CZ-gate. The x-part of the information flow vec- 
tor Xc(t) influences the measurement angles ±0 of the rota- 
tions for every 3-qubit gate according to Eq. ( 48 ) . As a side 



remark, the sign of the azimuthal angle 7r/2 of the measure- 
ment bases for the rotations also depends on the eigenvalues 
of the ancilla qubit according to Eq. ( 49 ) , which is not men- 
tioned in the table. 

We can easily generalize this example up the n-qubit gate 
A i2...(n-i) Z n Where the n - 1 logical qubits 1,2..., n - 1 
are the control qubits and the last one is the target qubit. 
For implementing this gate we need n — 2 work qubits, which 
are initialized in the state |+)®(™ -2 ). In Grover's search 
algorithm this gate plays a very important role, which we 
are discussing in the next section. 



IV. GROVER'S ALGORITHM WITHIN THE 
HQCM 

In this section, we illustrate the use of the HQCM in 
the context of a practical example: Grover's search algo- 
rithm (GA) [7]. Here, we shall see that the n-qubit gate 
^i2...(n-i)^n a i on g w ith the single-qubit operations is suffi- 
cient to perform GA. 

The best known quantum algorithm, for an unstructured 
database search, is GA. We have an unstructured database of 
total N items, out of which only one item matches with our 
query. Where a classical computer takes an average 0(N/2) 
steps, a QC with GA takes only &(VN) steps (Grover's it- 
erations) to find out the marked item. A "step" is a query 
of the oracle in the current context. 

A brief description of GA is given as follows: We can recog- 
nize each item of our database by a n-bit string. For simplic- 
ity, we choose N = 2™ = the total number of n-bit strings, 
and only one string out of N is marked. The task of the 
search problem is to recover the marked n-bit string in the 
end of computation. GA begins with the Hadamard opera- 
tion on each qubit, which were initially prepared in the state 
|0)® n , and that creates the superposition of the kets for all 
possible n-bit strings with equal amplitude. 

The next step is the implementation of Grover's iteration, 
which is a rotation in effect. It can be decomposed into two 
reflection operations. The first reflection operator is the or- 
acle O, which has the ability to recognize the solution of 
the search problem. Mathematically, O gives a conditional 
phase shift of tt to the matching string only. The second re- 
flection operator is the diffusion operator T>, which gives an 
inversion about the average. Like other quantum algorithms, 
GA is also probabilistic in nature. After 8(v / N) steps (iter- 
ations) the amplitude of the marked string becomes signif- 
icantly larger than those of the unmarked strings. Finally, 
we read the output by performing measurements on all the 
qubits. There have been many successful attempts at the im- 
plementation of GA, for N = 4 or N = 8, in different physical 
setups such as with NMR system [2TM23| . with cavity quan- 
tum electrodynamics (QED) [2lH2"5] . with optics [27J 125] . 
and with the MQCM [291 [30]. 

Now, let us take a look at the structure of these two re- 
flection operators. If only one out of N item matches with 
our query; then there exist N different O operations (one for 
each item). Mathematically, the oracle corresponds to the 
jth item can be written as 



The case of j = N, C N = 2"®" - 2|N)(N|, corresponds 
to |N) = and is nothing but the n-qubit gate 

7 Yi2...(n-i) _ a i reac jy discussed its implementation with 
the HQCM for n = 4 in Sec. |III C 3[ and its circuit diagram 
is shown in Fig. [4j Any other oracle can be derived by per- 
forming the gate(s) X on the relevant qubit (s) before and 
after performing the gate h}' 1 "^ n ~ 1 > Z n . For example, the 
oracle associated with the item |1) = |0)®" can be derived 
as 



X 



(52) 



So the gate A 12- "( n_1 ) Z n is used for implementing every or- 
acle. The information, "which of the oracle is executed by 
the black box?" , is hidden to us. In other words, we do not 
know on which qubit (s) the black box is implementing the 
X gate(s) along with A 12 -^ -1 ) Z n . 

The mathematical structure of the diffusion operator is 



V 



<8>n y®n 



-H W X 



12.(^1) 
(n-1) 



(53) 



0* =I® n -2\j){j\. 



(51) 



It can also be constructed by performing the Hadamard- and 
the A-gate on every logical qubit before and after performing 
the gate A 12 -^ n ~^ Z n . In summary, the gate A 12 -^" 1 )^" 
together with the Hadamard and the A gates is sufficient to 
realize GA in the framework of HQCM. 



V. SUMMARY AND OUTLOOK 

We have established the HQCM, which is a hybrid model 
of the MQCM and the UQCM, where at first, a big uni- 
tary gate under simulation is decomposed into a sequence of 
elementary gates. The elementary gates in the HQCM are 
an arbitrary single-qubit gate, the CZ gate and the multi- 
qubit rotation around the z-axis. Every single-qubit gate 
and the CZ-gate are realized by a unitary evolution. Ev- 
ery multi-qubit rotation is executed by preparing a respec- 
tive star graph state followed by a single measurement. The 
HQCM is a model where a portion of the quantum circuit 
is simulated by the unitary evolution, and the rest is by the 
measurements. 

The choice of elementary gates is governed by the exper- 
imental easiness in terms of resources and computational 
steps. The implementation of an arbitrary single-qubit gate 
with the unitary evolution is straightforward in comparison 
with its implementation with the MQCM. The CZ-gate in 
itself is the part of experimental setup for creating the graph 
states. The star graph states for the multi-qubit rotations 
can be realized in one shot, and a single measurement is 
enough for executing these rotations. 

The classical information processing in the HQCM is very 
simple in comparison with the MQCM. In the HQCM, only 
the 2n-component information flow vector and the propaga- 
tion matrices for the elementary gates are needed for the clas- 
sical information processing, and the total number of steps 
are taken by a classical computer in parallel for doing this 
is the total number of gates in a quantum circuit under sim- 
ulation. Furthermore, no preprocessing and no additional 
computational steps are required for classical information 
processing in the HQCM. 

We also have shown that how one can realize efficiently 
the multi-control gates (like A 12 -^ n ~ 1 '> Z n ) with the HQCM, 
which play a very important role in the implementation of 
GA. The gate A 12 -( n-1 >Z n together with the Hadamard and 



11 



the X gates is sufficient to realize GA in the framework of 
HQCM. 

One can carry on the investigation in the following direc- 
tions. In addition to the multi-qubit rotations in the set of el- 
ementary gates for the HQCM, one could include some more 
gates — which can be executed in one shot by the MQCM 
without adding further complications in the model — in the 
set of elementary gates. Further, implementation errors need 
to be considered for any practical realization. In order to es- 
tablish the fault-tolerance version of HQCM, one needs to 
design the elementary gates in a fault-tolerant manner but 
this is beyond the scope of the present paper. 



Appendix A: Single-qubit rotation Rzi'fi) 

In this appendix, we discuss the implementation of 
the single-qubit rotation around the z-axis Rzi'-P) = 
exp(~i(pZ/2) with the MQCM H2Q3]. A two-qubit graph 
state corresponding to the graph depicted in Fig. [5ji) is suf- 
ficient for accomplishing the job. The logical qubit 1 (rep- 
resented by the circle) is in a general input state |^i„(l)), 
and this is the single-qubit state on which we want to ap- 
ply R z (tp). In order to generate the required graph state, 
we prepare the qubit a (represented by the diamond) in the 
state (|0) o + (-1) K » |1)J /\/2. Then both the qubits are con- 
nected by the operation cz, which is represented by the bond 
in Fig. [5ji) and given by Eq. ([I]). The resulting graph state 



1 

71 



>a ® |^i„(l)} + (-1) K ° |l)a ® (^in(l)})] 

(Al) 

is ready for the computation. Here, the subscript 1+1 indi- 
cates that this graph state is made of two qubits, the logical 
qubit 1 and the ancilla qubit a. 

In order to generate the desired effect on the input state, 
we measure the qubit 1 in the basis 

{| t,l (tt/2, = {(10), + (-l)"*e-* (I),) 

(A2) 

and the value of mi is the result of the measurement. After 
the measurement, the output state (up to a global phase) 

|Vw(l)) = (AT 1 (Z)"' HR z (<p)\i,- m (l)) (A3) 

is obtained from the qubit a, and the qubit 1 gets projected 
either onto the state | f ( 7r /2, — f))i (if m i = 0) or onto 
the state | 4- (t/2, — cp))i (if mi = 1). The net effect on the 
input state is the required operation R z (<p) followed by the 
Hadamard gate 

H = (X + Z)/V2 

(represented by the boxes in Fig. [5jfi)) and the byprod- 
uct operator (A) 1 ™ 1 (Z) Ka (represented by the dotted-box in 

Fig- 6% 

Here, the direction of measurement lies in x,y plane of 
the Bloch sphere, and the input information is not only tele- 
ported from one lattice site to other but also gets processed 
by the measurement. As we know, every rotation in the 
Bloch sphere corresponds to a single-qubit operation (up 
to a global phase). Owing to the Eulcr decomposition of 
an arbitrary rotation R{a,P,^) = R z (j)R x (f3)R z (a), where 
Rx(fl) = exp (—i/3X/2), one can generate an arbitrary single- 
qubit operation with a chain of five qubits graph state with 
the measurement direction for each qubit (the angles a, /3, 7) 



(1) 



(ii) 



FIG. 5: (Color online) (i) The graph associated to the graph 
state |</>)(i + i) given by Eq. ( Al I. The (blue) circle, the bond, and 



the (black) diamond represent the logical qubit which is in the 
input state jf/>m(l)), the cz operations given by Eq. Q, and the 
ancilla qubit a, respectively, (ii) The quantum circuit illustrates 
the effect on the input state, when the qubit 'a' is measured in 
an appropriately chosen basis. 



in x,y plane of the Bloch sphere [THUS]- The realization of 
an arbitrary rotation by such a sequence of two z rotations 
sandwiching an x rotation illustrates the importance of the 
temporal ordering of the measurements in the MQCM. 



Appendix B: propagation matrices 

In this appendix, the propagation matrices for the R-, 
CNOT-, H-, and i? z (7r/2)-gate for the case of n logical qubits 
are given. The propagation matrix C is a 2n x 2n matrix of 



the form — — — , where C xx , C zx , C xz and C zz are 

V C xz C zz J 

nxn matrices with binary- valued entries |14) . One can gen- 
erate the propagation matrices for an arbitrary single-qubit 
rotation Ry> on the logical qubit c j' with 



c xx (R (j) ) 


kl 


C ZZ (R^) 


= 8 k i 

kl 




kl 


C XZ (RW) 


= 0; 

kl 



(Bl) 

Here, [C xx (R W )] , , stands for the entry in fcth row and Ith 
column of the matrix C xx corresponding to the R^ gate. 
The same notation apply for the propagation matrices de- 
fined below. 

The propagation matrix for the CNOT(a, b) gate (both the 
control qubit 'a' and the target qubit belong to the set of 
n logical qubits; a ^ b) is given by 

[C^CNOTta, b))] kl = S k i + S kb Si a , 
[C zz (cNOT(a, b))] u = S a + SkaSib, 

[C Z2; (CNOT( a , b))] kl = [C, z (CNOT(a, b))] kl = 0; (B2) 

for the Hadamard gate H^J on the logical qubit l f is given 
by 



(B3) 



and for the 7r/2-phase gate R^P {tt /2) on the logical qubit l f 
is given by 



C xx (HW) 


kl 


C ZZ (H^) 


= S k i © SkjSij 

kl 


~C ZX (HW) 


kl 


C XZ (H^)~ 


= SkjSij) 

kl 



C xx {R^{k/2)) 
C zx {R^{n/2)) 
C xz (R^(n/2)) 



ki 



C zz (R^(n/2)) 



= 5, 



kl 



kl ■ 



kl 



= SkiS, 



kl 



kjOly 



(B4) 



It is advantageous to deal with the information flow vector 
I together with the propagation matrices (Eqs. (B1|-(B4|) 



by a classical computer, than directly dealing with the cor- 
responding byproduct operator U b together with the prop- 
agation relations (Eqs. Jl4|, (fl5|), (fl8|, and (fl9|). 
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Appendix C: Efficient measurement scheme 

This appendix holds a discuss on an efficient measurement 
scheme of the MQCM where the temporal order of the mea- 
surements plays an important role [14j . On one hand in the 
UQCM, we cannot parallelize two gates of a sequence that 
do not commute. On the other hand in the MQCM, all the 
gates from the generating set of the Clifford group can be 
executed in a single time step irrespective of their positions 
in the circuit. In other words, the temporal order of the 
measurements in the MQCM is not pre-imposed by the tem- 
poral order of the gates. So, the most efficient scheme for 
the measurements does not necessarily follow the temporal 
order of the gates in a circuit under simulation. First, the 
spatial pattern of the measurement bases is assigned to the 
graph qubits according to the sequence of gates. Then the 
measurements are performed round by round according to 
the scheme which is given as follows. 

First, the graph Q is divided into disjoint subsets of qubits 
Qt , where index t stands for the round of measurements and 
< t < i max . Mathematically, Q = \J™g Q t and Q s nQ t = 
for all s 7^ t . The subset Qt is a collection of all those qubits 
which will be measured simultaneously in fth round. All the 
measurements in the A, Y and Z eigenbasis are put together 
in the very first round (Oth round), and there is no need for 
adjusting the measurement bases according to the previous 
measurement results for the qubits of Qq. In the first mea- 
surement round, the "redundant graph qubits" are removed 
by the Z-measurements, the "readout qubits" are measured 
in the Z eigenbasis, and the Clifford part of the circuit is ex- 
ecuted by the X-, Y-measurements [T3]. In the MQCM, the 
"readout qubits" which play the role of "output register" are 
not the last ones to be measured, they are among the first 
ones. 

The observables for all subsequent measurement rounds 
are of the form cos(ip)X ± sm(<p)Y with ip ^ {0, ±f }, where 
the measurement bases for the qubits are decided by the 
measurement outcomes from the previous rounds. All those 
qubits whose measurement bases depend on the outcomes 
from the first measurement round belong to the next subset 
Q\. Similarly, the measurement outcomes from the subset 
<2i together with Q decide the measurement bases for the 
qubits in Q 2 , and so on. These subsets are measured one by 
one up to the final measurement round i max . One can think 
of the total number of measurement rounds (i max + 1) as the 
logical depth (temporal complexity) for the MQCM. 

Parallel to the measurement rounds, the classical data- 



processing parts are taken care of by a classical computer. 
After preparing the graph state and just before starting the 
measurements, the information vector is initialized to 2^? t 0M - 
^in?™ depends on the eigenvalues {k} of the graph and some 
particular gates (like CNOT, R z (n/2)) which appear in a 
quantum circuit under simulation |14j . After executing the 
first measurement round on the set Qq, 2^?t° M gets updated 
to Z(0) through the measurement results. 1(0) then deter- 
mines the measurement bases for the qubits of Q%. Similarly, 
the measurement outcomes from round t update the infor- 
mation flow vector from X(t — 1) to I(t). The corresponding 
byproduct operator is given by 

n 

u B {t) = Y[( xW )*' (t \z U) ) Zi( - t) ■ (CI) 
J=l 

Then, I(t) = I(xj(t), Zj(t)) sets the measurement bases for 
the (i + l)th round. After the final measurement round t max , 
the a;-part of the information flow vector Xr(i m ax) enables us 
to interpret the result of the computation. 

In this measurement scheme, the following technical points 
are worth emphasizing, which are discussed in Ref. |14j : (1) 
In order to construct the subsets of graph qubits {Qt}, a 
classical computer needs the forward cones for all the graph 
qubits. The forward cones decide a strict partial ordering 
among the qubits, and the sets {Qt} are constructed accord- 
ingly. (2) In order to account for the influence of the mea- 
surement outcomes and the set of eigenvalues {k} on I(t), 
a classical computer needs the byproduct images for all the 
graph qubits. (3) {k}, the byproduct images, and I(t) are 
required for setting the measurement bases for the as-yet un- 
measured qubits. A classical computer uses the symplectic 
scalar product for doing this. 
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